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Q , ABSTRACT 

Multifrequency observations with the Giant Metrewave Radio Telescope (GMRT) and 
the Very Large Array (VLA) are used to determine the spectral breaks in consecutive 
strips along the lobes of a sample of selected giant radio sources (GRSs) in order to 
estimate their spectral ages. The maximum spectral ages estimated for the detected 
\ radio emission in the lobes of our sources range from ~6 to 36 Myr with a median 

<*S |. value of ^20 Myr using the classical equipartition fields. Using the magnetic field 

' ' estimates from the Beck & Krause formalism the spectral ages range from ^5 to 38 

Myr with a median value of ^22 Myr. These ages are significantly older than smaller 
sources. In all but one source (J1313+6937) the spectral age gradually increases with 
^3 I distance from the hotspot regions, confirming that acceleration of the particles mainly 

occurs in the hotspots. Most of the GRSs do not exhibit zero spectral ages in the 
hotspots, as is the case in earlier studies of smaller sources. This is likely to be largely 
due to contamination by more extended emission due to relatively modest resolutions. 
The injection spectral indices range from ^0.55 to 0.88 with a median value of ^0.6. 
We discuss these values in the light of theoretical expectations, and show that the 
injection spectral index appears to be correlated with luminosity and/or redshift as 
■ well as with linear size. 

£NJ ■ Key words: galaxies: active - galaxies: jets - galaxies: nuclei - quasars: general - 

radio continuum: galaxies 
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1 INTRODUCTION has been seen in many sources and used to estimate the ra- 

diative ages and expansion velocities in several samples of 
♦H The radio continuum spectra m different parts or an ex- „ , „„„ / „ „ , ,„ or A1 

rrt , , . . . - . , powerful 3CR sources (e.g. Myers & Spangler 1985; Alexan- 

. . 1 tended radio source contain important information about , „ T , , „„- T 1 1, 1 r., 1 -1^0^ ^ -n- 

. , , ... . . der & Leahy 1987; Leahy, Muxlow & Stephens 1989; Canlh 

the various energy losses and gams 01 the radiating particles , ,„_., T . , n-i i™r,\ • , ri 

, r r , , , , , et al. 1991; Liu, Pooley & Riley 1992), m samples 01 low- 

durmg the lifetime of the source. According to dynamical , , ,. , ,. , . , r ., . 

1 r nnTT n- o n n luminosity and medium- luminosity radio galaxies (e.g. Klein 

models lor FRII-type (Fanaroft & Riley 1974) radio sources, , lrl __ . „ , r 

, . . . et al. 1995; Parma et al. 1999), as well as a sample ot com- 

the energy emitted by the nucleus squirts out m the iorm , „, . , -ir,r,r,\ tt 

... ,. , .. .. . . . pact steep-spectrum sources (Murgia et al. 1999). However, 

or narrow, colhmated jets channeling their way through the 



external environment. These jets dissipate their energy at 



the observed spectral break and steepening of the spectrum 

beyond this break need not be entirely due to radiative en- 
their leading edges giving rise to intense regions of emission , . ., , , r , , , , ,. n , , , 

, a 1 f , er gy fosses. A possible role of the focaf magnetic fiefds, de- 

calfed hotspots . As the lets advance outwards, the refativis- , rjl , ln rjl , , . , jn ,.„, nj . 

tails of the backftow ot the fobe material, or the difhcufties 
tic particles now out from the hotspots to form the extended ,. ,. . , „ jn .... 

, m disentangling various energy losses ot the radiating par- 

lobes of radio emission. Assuming that there is no significant ,-11 1 ■ , , . - , r 1 

ticfes have been pointed out m a number ot papers (e.g. 

reacceleration within these lobes and no significant mixing t,t--, ^ 1 • 1 1™,-, n 1 • 1 is l. n, c a 

, , , , , Wuta & Gopal-Krishna 1990; Rudmck, Katz-Stone & An- 

of particles, there should be a spectral gradient across the , , , , „ „ „ , 

, ' . , , , 1 derson 1994; Eilek & Arendt 1996; Jones, Ryu & Engel 1999; 

radio source, m the sense that the spectrum should steepen 1 n c n t mm \ 

... . ,. r . _,. . ... Blundell & Rawlmgs 2000). 

with increasing distance from the hotspot. this prediction 

Nevertheless, giant-sized radio sources (hereinafter re- 
ferred to as GRSs) are suitable for the classical spectral- 
* E-mail: jamrozy@oa.uj.edu.pl (MJ); chiranjib@iucaa.emct. in ageing analysis due to their large angular extent which can 

(CK); machalsk@oa.uj.edu.pl (JM); djs@ncra.tifr. res. in (DJS) be covered by a significant number of resolution elements. In 
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Figure 1. Upper panel: an example of how the lobes of analysed sources are cut into strips within which the spectral age of radiating 
particles is determined. Lower panels: the spectra fitted to the flux-density data in each of the strips. 



Paper I (Konar et al. 2007) extensive low-frequency observa- 
tions made with the GMRT of a sample of ten selected GRSs 
have been described and some of their properties analysed 
and discussed. In this paper we analyse the spectral ages 
of these radio galaxies except for J0720+2837 for which we 
have data at only two frequencies and J1343+3758 whose 
spectral ageing analysis has been presented by Jamrozy et 
al. (2005) . Our analysis is similar to the spectral ageing anal- 
ysis for the GRS J1343+3758 (Jamrozy et al. 2005) and the 
double-double radio galaxy J1453+3308 (Konar et al. 2006). 
In Section 2 we discuss briefly the method of analysis, while 
in Section 3 we present our results. Discussion of the results 
and an examination of the correlations of the injection spec- 
tral index, ainj, with luminosity, redshift and linear size are 
presented in Section 4. The conclusions are summarized in 
Section 5. 



2 SPECTRAL AGEING ANALYSIS 

2.1 Spectral steepening and spectral age 

The spectral age, i.e. the time which has elapsed since the 
particles were last accelerated in the hotspot region, can be 
derived from the steepening in the radio spectrum due to 



energy losses caused by synchrotron and inverse- Compton 
processes. The observed spectra have been fitted using the 
Jaffe & Perola (1973, JP) and the Kardashev-Pacholczyk 
(Kardashev 1962; Pacholczyk 1970; KP) models using the 
SYNAGE package (Murgia 1996). We found no significant dif- 
ference between these two models over the frequency range 
of our observations and present only the results obtained 
from the JP model in this paper. Also, while fitting the the- 
oretical spectra to the observed ones, we find the spectra of 
the different strips in the lobes are better fitted with the JP 
model than with the continuous injection (Kardashev 1962; 
CI) model, though sometimes the CI model fits better in the 
area of a prominent hotspot. 

Thus, under assumptions that (i) the magnetic field 
strength in a given lobe is constant throughout the energy- 
loss process, (ii) the particles injected into the lobe have a 
constant power-law energy spectrum with an index 7, and 
(iii) the time-scale of isotropization of the pitch angles of 
the particles is short compared with their radiative lifetime, 
the spectral age, r spec , is given by 

r spcc = 50.3 K,(l + AY 1 ' 2 [Myr], (1) 
where f?ic=0.318(l+z) 2 is the magnetic field strength equiv- 
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Figure 2. Spectral-index map, flux-density profiles, and the spec- 
tral age distribution for J0912+3510. The spectral ages have been 
estimated using magnetic field values determined using the Beck 
& Krausc (2005) formalism (filled circles and with error bars) and 
the classical (e.g. Miley 1980; referred to as classical-1 in Paper 
I) formalism (open circles without error bars). 



alent to the cosmic microwave background radiation. Here 
B, the magnetic field strength of the lobes, and Bic are ex- 
pressed in units of nT, v br is the spectral break frequency 
in GHz above which the radio spectrum steepens from the 
initial power-law spectrum given by a in j=(7 — l)/2. It may 
be noted that Alexander & Leahy (1987) and Alexander 
(1987) have shown that the effects of expansion losses may 
be neglected. 



2.2 Determination of a^j and v br values 

In order to determine a value of Qinj, we fit the JP model 
to the flux densities of the entire lobes (given in Paper I) 
treating ai n j as a free parameter. We find that the uncer- 
tainties of the ai n j values are sometimes large due to the 
low flux densities of the lobes. We have used the values of 
Qinj which give the best fits to the spectra of the lobes, and 
note that there is usually no evidence for significantly dif- 



Table 1. Break frequency, magnetic field strength, and spec- 
tral age of particles in consecutive strips through the lobes for 
J0912+3510 



Strip Dist. 

(kpc) 



(GHz) 



S eq (rev) 
(nT) 



' spec 

(Myr) 



PI 
P2 
P3 
P4 
P5 
P6 



P7 

P8 

P9 

P10 

Pll 

P12 



543 
466 
389 
313 
236 
161 



498 
575 
652 
729 
806 
883 



N-lobe 

+32.9 

7.7 
>+85 
11.8 
>+90 
12.5 
+48.0 
-4.9 
>+70 
3.5 

4.9>+f 



12.6 
17.2 
18.6 
8.8 
8.6 



6.7 
8.1 
14.9 
21.6 
39.6 
22.6 



S-lobe 

>+90 

1.9 
+38.6 

4.3 
+ 13.2 
10.7 
+35.2 
15.9 
+96.7 
31.3 
+23.0 
16.7 



0.07 
0.02 
0.35 
0.08 
0.07 
0.38 



9.93 
0.00 
2.05 
1.06 
2.52 
6.92 



ct inj = 0.560 
0.39±.032 
0.31±.028 
0.31±.027 
0.29±.027 
0.26±.025 
0.22±.022 

a inj = 0.628 
0.28±.031 
0.31±.034 
0.34±.037 
0.40±.042 
0.41±.042 
0.45±.046 



20.1 
17.9 
17.2 
25.1 
25.2 
32.6 



+9.8 

12.1 
+ 14.0 
17.0 
+ 12.9 
16.3 
+ 12.5 
-15.3 
+7.8 

23.8 
+ 14.5 
32.0 



7+5.3 
-28.6 
-12.1 
-15.3 
-,+5.4 

-16.9 
.,+6.2 
-14.4 
n _ 2 +5.3 



28. V_ 

26. i; 1 
19.0j 
15. 3^ 



14.3 



10.3 
+4.6 
13.8 



ferent values of this parameter for the oppositely directed 
lobes. For sources with reliable flux densities at lower fre- 
quencies than our observations but with no information on 
the flux densities of the individual lobes at these frequen- 
cies, we have determined ai n j from the integrated spectra 
of the sources. For these cases we have used the same value 
of Qinj for both the lobes. Having estimated the ai n j values, 
the total-intensity maps of the GRSs published in Paper I 
are convolved to a common angular resolution, which corre- 
sponds to the lowest resolution of our images. Each lobe is 
then split into a number of strips, separated approximately 
by the resolution element along the axis of the source in such 
a way that the extreme strips are centred at the peaks of 
brightness on the convolved maps. In each source a region 
around the radio core is excluded. As an example, the divi- 
sion of the source J1604+3731 into different strips is shown 
in Fig. 1. Using the SYNAGE software we determine the best 
fit to the spectrum in each strip over the entire observed 
frequency range using the JP model, and derive the relevant 
value of i^br- The fitted spectra in the strips covering the 
SE-lobe and the NW-lobe of J1604+3731 are plotted in the 
lower panels of Fig. 1. 



2.3 Magnetic-field strength determination 

In order to estimate the spectral ages in different parts of 
the lobes, we have to estimate the magnetic-field strength in 
the consecutive strips. Following Konar et al. (2006, 2007), 
the values of the equipartition energy density and the cor- 
responding magnetic field, B oq , are calculated using the re- 
vised formalism proposed by Beck & Krause (2005, here- 
inafter referred to as BK). This has been described in some 
detail in Paper I, and the values are denoted by B e q(rev) 
in this paper. The revised field strength within each of 
the consecutive strips is calculated assuming a filling fac- 
tor of unity. Since the BK formalism yields magnetic field 
strengths which are larger by a factor of ~3 compared with 
the classical formalism (e.g. Miley 1980; hereinafter refer- 
erred to as classical-1, as in Paper I), we also estimate the 
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Table 2. As in Table 1 but for J0927+3510 
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Figure 3. As in Fig. 2, but for J0927+3510 



spectral ages using the classical- 1 magnetic field strengths. 
This shows the range of ages possible for different estimates 
of the magnetic field strength. Parameters derived using the 
latter field strenghts are denoted by 'class' in this paper. It is 
relevant to note that in most of the previously published pa- 
pers on spectral ageing analysis of smaller FRII-type radio 
sources the spectral ages were determined using the classical 
minimum-energy formalism. 



3 THE OBSERVATIONAL RESULTS 

The results of our spectral ageing analysis for the eight GRSs 
selected from Paper I are presented in Figures 2 to 9 while 
some of the parameters are tabulated in Tables 1 to 8. In 
Figs. 2 to 9, (i) the spectral index map using the lowest 
frequency GMRT data and the VLA 4.86 GHz data, (ii) the 
intensity profiles of the lobes measured along the source axis 
at the different frequencies, and (iii) the resulting spectral 
age as a function of the distance from the radio core, are 
shown in the top, central and bottom panels respectively. 
Some of the numerical values are given in Tables 1 to 8 
which are arranged as follows. Column 1: identification of 



Strip Dist. 

(kpc) 



(GHz) 



^red 



Beq(rev) 
(nT) 



Tspcc 

(Myr) 



PI 
P2 
P3 
P4 
P5 



P6 
P7 
P8 
P9 



1125 
934 
742 
550 
217 



262 
665 
857 
1049 



SE-lobe 

20.6>+f 

rj + 85.0 

-15.9 
-, + (,1.6 

-5.7 
n +22.2 

-2.5 
I >+90 
-4.2 



22.2Z 

10.0I 1 

4.0^ 
8.1 ; 



NW-lobe 

8.4+f n 4 



35.2 
188.8 



2.0 

>+90 
11.9 

>+85 
134. 



0.36 
0.02 
0.05 
0.01 
0.96 



0.14 
2.89 
9.73 
4.19 



Q inj = 0.700 
0.47±.140 
0.28±.099 
0.24±.077 
0.20+.066 
0.21+.069 



0.25+.116 
0.34±.151 
0.37+.166 
0.40+.175 



7+2.5 

-5.5 

3 +5.4 
-6.5 

11.8+ 4 ' 4 



7.7Z 
6.9"* 



14.6; 
10.4" 



8.1 

8.9 
11.2 
6.0 
10.0 



14.5 

11.7 
5.8 
2.5 



+5.9 
6.0 
+2.1 
-8.2 
1.4 
5.5 
+2.2 
2.4 



J1 155+4029 
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Figure 4. As in Fig. 2, but for J1155+4029 
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Table 3. As in Table 1 but for J1155+4029 



spectral index 605 MHz - 4873 MHz 



Strip 


Dist. 
(kpc) 


(GHz) 


Y 2 


Boq(rev) 
(nT) 


Tspcc 

(Myr) 


PI 
P2 


925 
812 


SW-lobe 

1 93 1 + 85 -° 

40 7>+ 90 
4U - ' -32.9 


0.17 
0.25 


a inj =0.838 
0.84 ± .275 
0.68 ± .231 


9 9 +2.1 
5 2+ 4 ' 2 


P3 
P4 
P5 


226 
346 
459 


NE-lobe 

6.3t™ 

n 7 >+90 

81 S>+ 90 
81 - S -34.6 


0.06 
1.17 
1.26 


a inj =0.876 
0.93 ± .115 
1.23 ± .145 
1.62 ± .184 


11 1+ 3 ' 6 
q 1+0.9 

1 s+0-6 
1 - 8 -1.5 


Table 4. As in 


Table 1 but for J1313+6937 




Strip 


Dist. 
(kpc) 


(GHz) 


V 2 
A red 


Beq(rev) 
(nT) 


7"spec 

(Myr) 



PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

P9 

P10 

Pll 



470 
430 
390 
349 
309 
268 
228 
188 
147 
107 
67 



SE-lobe 

8.0±1;« 



"in 



: 0.610 



>+80 
1.5 

14.5>+»° 

l>+80 
-10.5 
7+40.0 
-28.9 
3 +34.5 
-25.3 
9 o c + 17.7 
23-»_ 186 
!>+80 
-17.1 

12.5+^ 

io.iti.S 

,+3.2 
-7.3 



7.5 



19.4; 
34.7^ 
30.9"* 



22.9; 



11.3"* 



0.73 

I. 82 

II. 1 
0.07 
4.38 
3.29 
3.87 
0.42 
0.16 
6.34 
4.25 



0.47±.015 
0.44±.014 
0.46±.015 
0.47±.015 
0.48±.015 
0.50±.016 
0.47±.015 
0.40±.014 
0.38±.013 
0.46±.015 
0.55±.017 



31 4+ 5 - 8 
01 -^-15.4 

34 1+ 4 ' 3 
23 9+ 3 - 2 
20 3+ 9 ' 7 
14 9+ 4 ' 9 

i0 -°-15.0 

18 4+ 4 ' 8 
lo -^-18.0 
+ 12.3 
20.3 
+ 10.5 
14.0 
+3.4 
9.8 
+3.3 
16.2 



20.5 
28.7 
28.6 
23.7 







NW-lobe 




« inj = 0.610 




P14 


31 


8-131 


0.38 


0.54±.017 


28.0 


P15 


71 


9-Otjjl 


1.22 


0.51±.016 


28.0 


P16 


111 


„ q>+80 


0.27 


0.48±.015 


26.2 


P17 


152 


15.2>+f 


1.03 


0.50±.016 


21.9 


P18 


192 


22.5 >+ 3 5 


0.99 


0.55±.017 


16.7 


P19 


230 


24 Q>+85 
z ^- y -15.9 


0.001 


0.59±.018 


15.1 


P20 


272 


28.4+ 42 ;« 


1.55 


0.65±.019 


12.7 


P21 


313 




0.96 


0.71±.020 


12.6 



+4.9 
-6.0 
+4.9 
-15.3 
+6.8 
-20.6 
+5.0 
-19.1 
+5.1 
-15.9 
+9.9 
-13.4 
+4.9 
-12.5 
+5.2 
-12.4 



the strip, column 2: projected distance of the centre of the 
strip from the radio core in units of kpc, column 3: break 
frequency of the spectrum of the strip according to the JP 
model in units of GHz, column 4: reduced \ 2 value of the 
fit, column 5: revised magnetic-field strength in units of nT, 
and column 6: spectral age of particles in the strip. We have 
estimated the physical parameters using the cosmological 
model where H D =71 km s _1 Mpc -1 , f2 m =0.27, f2„ ac =0.73 
(Spergel et al. 2003). 



3.1 Notes on individual GRSs 

J0912+3510: there is some uncertainty as to which galaxy is 
associated with the radio source (cf. Machalski et al. 2006, 
hereinafter referred to as MJZK). 

J0927+3510: as mentioned in Paper I, the strips P5 
and P6 in Table 2 overlap the region of emission which was 
suspected to be an inner double (cf. MJZK). However, the 
very steep spectra of about 1.17 and 1.35 in the strips P5 
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Figure 5. As in Fig. 2, but for J1313+6937 



and P6 respectively suggest that this is likely to be old lobe 
emission. The resulting spectral ages are not significantly 
different from the expected trend of age with distance from 
the hotspot, with the age of P5 which is weaker and on the 
south-eastern side being only marginally lower. The some- 
what higher separation speed of ~0.7 and 0.4c for the east- 
ern and western lobe material respectively (Fig. 3) suggests 
that the photometric-redshift estimate of 0.55 for this GRS 
(MJZK) may be overestimated. 

J1155+4029: although we have tentatively identified it 
as a galaxy (MJZK; Paper I), an optical spectrum of the 
faint host would help resolve the nature of the optical galaxy. 

J1313+6937, DA 340: the spectral-index map and the 
intensity profiles along the source axis suggest that the syn- 
chrotron radiation observed in separate strips of this GRS 
may be related to a mixture of emitting particles which were 
injected or accelerated at different time epochs. The SE lobe 
of this GRS is the only one, among the other lobes studied 
in this Paper, in which the synchrotron ages determined do 
not more or less gradually increase with distance from the 
hotspot area to the core. The spectral ageing analysis of 
the entire lobes of this galaxy was published by Schoenmak- 
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Figure 6. As in Fig. 2, but for J1604+3438 



Figure 7. As in Fig. 2, but for J1604+3731 



Table 5. As in Table 1 but for J1604+3438 



Table 6. As in Table 1 but for J1604+3731 



Strip 


Dist. 
(kpc) 


(GHz) 




Beq(rev) 
(nT) 


7~spcc 

(Myr) 






E-lobe 




a inj = 0.554 




PI 
P2 
P3 
P4 


441 
376 
312 
248 


26.9+«g-; 
33.6>«° 

»*±\\ 
4 3+1-0 


0.82 
0.67 
0.98 
0.16 


0.64±.021 
0.52±.015 
0.50+..014 
0.53+..016 


10 1+ 4 ' 8 

10 2+ 4 - 7 

19 4+ 101 
la - -2.9 

28.2+3 \ 


P5 
P6 


183 
119 


2 5+ 12 

Z -°-0.3 
n K+4.9 
Z -°-0.5 


0.26 
1.87 


0.51+..014 
0.57±.017 


37.5±*» 
35-8+t ! 







W-lobe 






a in j = 0.554 






P7 


87 


6 >+9 ° 
D ' u -0 .7 


1 


39 


0.62±.020 


21 


9 


P8 


198 


3 6+ 4 1 





93 


0.47+..013 


32 


5 


P9 


262 


6.0+i^ 





71 


0.46+..012 


25 


8 


P10 


327 


11 4+ 3 ' 3 


1 


14 


0.47+..013 


18 


3 


Pll 


391 


20.2+25j 





39 


0.50+..014 


13 


4 



+ 1.9 
-21.0 
+2.6 
-10.3 
+3.9 
-4.9 
+ 11.6 
-2.2 
+4.5 
-9.0 



Strip Dist. 

(kpc) 



PI 
P2 
P3 
P4 



(GHz) 



(nT) 



S-lobe 



647 
541 
435 
329 



31.7 



>+80 



20.0 
8.6 



18.4 
>+90 
10.6 
+41.8 
-2.0 
4 2 +6.9 
^-0.6 
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Figure 8. As in Fig. 2, but for J1702+4217 
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Figure 9. As in Fig. 2, but for J2312+1845/3C457 



Table 7. As in Table 1 but for J1702+4217 
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Table 8. As in Table 1 but for J2312+1845 
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ers et al. (2000) where the radio spectrum of each of the 
lobes were determined within the frequency range of 325 
MHz-10.5 GHz. Though their ai n j values found fitting the 
JP model, i.e. 0.71±0.02 and 0.75±0.02 for the SE and NW 
lobes, respectively, differ from our values, their average syn- 
chrotron ages of these lobes of 28±1 Myr and 27±1 Myr are 
consistent with our results (cf. Table 4). 

J1604+3438: the inner part of the structure comprising 
the fiat-spectrum radio core and the two nearby peaks of 
emission seen most clearly in the image published by MJZK 
led them to classify the source as a possible double-double 
radio galaxy (DDRG). However as noted in Paper I spectral 
indices of these features are steep and one needs to image 
them with higher resolution in both total intensity and linear 
polarisation to investigate the DDRG nature of the source. 
It is relevent to note that the spectral age of the western 
feature lies significantly below the age— distance relationship 
for the western lobe. 

J1604+3731: this is the highest-redshift radio galaxy 
in the sample. The relatively younger age seen in the strip 
closest to the core in the northern lobe is coincident with a 
knot or peak of emission (Konar et al. 2004; Paper I). 

J1702+4217: although there are peaks of emission vis- 
ible in the lobe emission, especially at the lower frequency, 
the spectral age increases smoothly with distance from the 
hotspots. 

J2312+1845, 3C457: the only source in the sample 
showing undoubtedly a curvature in the age— distance plot. 
This could be due to deceleration of the radial expansion 
velocity of the lobe material. The data in Table 8 give an 
expansion velocity of about 0.11c±0.05c in the vicinity of 
the leading heads which slows to about 0.03c±0.015c in the 
regions where the oldest particles are detected. Such curva- 
tures may also be produced by re-acceleration of particles in 
different regions of the lobes. 

3.2 Separation speeds and spectral ages 

The plots in the bottom panels of Figs. 2-9 show that, in 
general, the spectral age of the emitting particles increases 
systematically with distance from the leading heads of the 
lobes, i.e the hotspot regions. These plots have been shown 
using the magnetic field strengths estimated using the BK 
formalism (full circles) as well as the classical-1 (Milcy 1980; 
see Paper I) minimum-energy one (open circles) . These two 
field strengths differ by a factor of ~3 (cf. Paper I). The 
numerical values we list here have been estimated using the 
BK formalism. For most of these lobes we calculate the linear 
regression of the ages (r spec ) on the corresponding distances 
(d), which is shown as a dashed line in Figs. 2-9. The slopes 
of these lines give a characteristic speed, « sep , which is an 
indication of the average speed of the lobe material relative 
to the hotspots. A clear exception is J2312+1845 (3C457) 
for which the speed, v scp , evidently decelerates with distance 
from hotspot regions. This source is much stronger than the 
remaining GRSs, therefore the errors of the spectral break 
frequency, v br , and the corresponding spectral age in each 
of the strips are much smaller than those in our remaining 
sources. In other cases the linear regression fit is the most 
conservative, though there is no physically justified reason 
for the material separation speed to be constant. 

Our estimates of the average speed and a 'characteristic' 
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Figure 10. The spectral age in Myr as a function of the largest 
linear size in kpc 
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spectral age for the lobes in the regions of the observed emis- 
sion from the lobes as well as the values when extrapolated 
to the positions of the cores are presented in Table 9. These 
values have been presented for all the sources except for the 
SE-lobe of J1313+6937 which has a complex age— distance 
plot. In Table 9 the values have been presented for the esti- 
mates using the magnetic field values from both the BK and 
the classical-1 formalisms. The Table is arranged as follows. 
Column 1: source name; column 2: the optical identification; 
column 3: redshift; column 4: largest linear size / in kpc; col- 
umn 5: log of the luminosity in units of W Hz -1 ; column 6: 
lobe identification; column 7: spectral age using the BK for- 
malism, denoted by 'rev', till the most distant strip from 
the hotspot; column 8: spectral age when the age-distance 
relationship is extrapolated to the region of core; columns 9 
and 10: same as columns 7 and 8 but using the classical-1 
estimates of the magnetic field, denoted by 'class'; columns 
11 and 12: estimates of the separation velocity for the BK 
and classical-1 approaches. 
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Table 9. Some of the physical parameters of the sources 
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25.44 


N 


28.7 


49.9 


23.0 


34.2 


0.09 


0.18 
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32.6 


35.1 


22.2 


23.2 


0.07 


0.12 


10927+351 n 
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0.55* 


2206 


26.02 


SE 


14.6 


21.4 


8.3 


10.4 


0.24 


0.70 












NW 


14.5 


19.6 


8.7 


12.2 


0.21 


0.37 


T1 1 tt-i-dftVQ 
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0.53* 


1437 


26.55 
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0.12 


0.12 
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0.08 


j loioTuyo / 
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NE 
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0.13 
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52.5 
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0.12 


0.09 
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0.249 
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S 


58.4 


65.0 


47.9 


61.8 


0.04 


0.05 












N 


46.6 


54.9 


38.8 


53.1 


0.04 


0.05 


J1604+3438 


G 


0.2817 
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25.53 


E 


56.3 


56.3 


34.7 


50.8 


0.03 


0.03 
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52.1 


28.7 


44.5 


0.03 


0.04 


J1604+3731 


G 


0.814 


1346 


26.60 
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15.8 


15.8 


7.2 


11.3 


0.16 


0.25 
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0.24 
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1160 
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16.3 
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0.12 


0.16 
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0.18 


J2312+1845 


G 


0.427 


1056 


27.11 


NE 


34.8 


34.8 


21.9 


36.5 


0.08 3 


0.07 3 












SW 


49.2 
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Notes: * denotes an estimated redshift; (1) calculated from the data in Jamrozy et al. (2005); (2) calculated from the data in Konar et 

al. (2006); (3) speed estimated using the two extreme strips for each lobe. 



4 DISCUSSION 
4.1 Spectral ages 

The maximum spectral ages estimated for the detected ra- 
dio emission in the lobes of our sources range from ~6 to 
36 Myr with a median value of ~20 Myr using the classical 
cquipartition fields. Using the magnetic field estimates from 
the BK formalism the spectral ages range from ~5 to 38 Myr 
with a median value of ~22 Myr. These ages are significantly 
older than smaller sources. The median linear size of these 
large sources is ~1300 kpc. For comparison, we first con- 
sider the sample of smaller sources studied by Leahy et al. 
(1989) which has been observed over a large frequency range. 
They estimated the spectral ages of a sample of 16 3CR ra- 
dio sources from low-frequency observations with the Multi- 
Element Radio Linked Interferometer Network (MERLIN) 
at 151 MHz and high-frequency observations at 1500 MHz 
with the VLA. We have re-estimated their values using the 
cosmological parameters used in this paper. The ages of the 
lobes of these sources range from ~2.5 to 26 Myr with a me- 
dian value of ~8 Myr, significantly smaller than the value 
for the GRSs. The median value of the linear size for this 
sample of sources is 342 kpc. 

Considering the sample of 14 high- luminosity double 
radio sources studied by Liu et al. (1992), we have again re- 
estimated the spectral ages using the cosmological parame- 
ters used in this paper, and find that the values range from 
~0.3 to 5.3 Myr with a median value of ~1.7 Myr. These 
sources have even smaller spectral ages and their median 
linear size is 103 kpc. Considering these different samples, 
there is a trend for the spectral ages to increase with linear 
size (Fig. 10) as has been noted earlier in the literature (e.g. 



Parma et al. 1999; Murgia et al. 1999; Murgia 2003). It is 
also interesting to note that the relative speeds of the lobe 
material has been estimated to be in the range of ~0.03 to 
0.2c with a median value of ~0.1c for the Liu et al. sam- 
ple. In comparison, the corresponding speeds for our sample 
of GRSs is about 0.03 to 0.25 (excluding J0927+3510 for 
which the photometric redshift value may be overestimated, 
cf. Sect. 3.1), with a median value of ~0.09c, similar to that 
of the the Liu et al. sample. 

The above trend for spectral age to increase with size is 
broadly consistent with the expectations of dynamical mod- 
els of the propagation of jets in an external medium (e.g. 
Falle 1991; Kaiser & Alexander 1997; Jeyakumar et al. 2005 
and references therein) in which the linear size of the source 
is a function of its (dynamical) age. Nevertheless, while in- 
terpreting these numbers caveats related to the evolution 
of the local magnetic field in the lobes need to be borne in 
mind (e.g. Rudnick, Katz-Stone & Anderson 1994; Jones, 
Ryu & Engel 1999; Blundell & Rawlings 2000). Also, while 
Kaiser (2000) have suggested that spectral and dynamical 
ages are comparable if bulk backflow and both radiative and 
adiabatic losses are taken into account in a self-consistent 
manner, Blundell & Rawlings (2000) suggest that this may 
be so only in the young sources with ages much less than 
10 Myr. In the study of the FRII type giant radio galaxy, 
J1343+3758, Jamrozy et al. (2005) find the dynamical age 
to be approximately 4 times the maximum synchrotron age 
of the emitting particles. 

In the plots of spectral age against distance from the 
hotspots, we note that none of the sources show zero age at 
the hotspots. In fact the ages in the hotspots usually range 
from a few to ~20 Myr, except in the case of the south- 
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eastern lobe of J1313+6937 where the age in the outer ex- 
tremity is ~30 Myr. This tendency for the hotspots to ex- 
hibit a non-zero age has been noted earlier by a number of 
authors (e.g. Liu et al. 1992). In the Liu et al. sample, ages 
in the hotspots are typically less than ~1 Myr while in the 
case of Cygnus A, Carilli et al. (1991) find the ages of the 
hotspots to be 0.1 and 0.2 Myr for the north-western and 
south-eastern lobes respectively. The typical linear resolu- 
tion of the sources observed by Liu et al. is ~10 kpc for most 
sources while that of Cygnus A is ~5 kpc. In comparison, 
the linear resolution of our GRSs ranges from ~30 to 110 
kpc with a median value of ~80 kpc. For comparison, the 
sizes of hotspots range from about a few parsec to ~10 kpc 
increasing with source size till the largest linear size is ~100 
kpc and then flattens for larger source sizes (e.g. Jeyakumar 
& Saikia 2000). Possible explanations for the non-zero ages 
include (i) the possibility that the hotspots have been inac- 
tive for 10 5 to 10 6 yr, (ii) contamination by more extended 
emission and (iii) higher magnetic fields in the hotspots (cf. 
Liu et al.). The coarse resolution of our observations suggest 
that contamination by extended emission is the likely cause 
for the non-zero ages in the hotspots of our GRSs. 

Another effect, visible in the bottom panels of Figs. 2-9 
and related to the magnetic fields, is the age difference, A = 
T S pec(rev) — r spoc (class), which is either positive or negative. 
This is a simple consequence of Eq. (1). Given Ubi and z 
imply the maximum Tspec 

for B=Bic/V3. Substitution of 
the ratio r = B /B\c into Eq. 1 gives 

Tspec = 50 - 3 Y^ #- 3/2 {Ml + Z)V 1/2 = 

= 50.3fc(r){0.318(l + z) 2 y i/2 {v hl {l + z)}" 1/2 . 

The function k(r) has a maximum at r*=3~ 1 / 2 , thus the 
maximum of r spec can be written as 

r s ( p m c ? x) = 159.8 v£'\l + z)- 7 ' a 

The function k(r) as a function of r is plotted in Fig. 11. 

In our calculations the values of the magnetic field 
calculated with the revised formula, B oq (rev), are always 
greater than the values of B eq (class) calculated with the 
classical formula of Miley (1980) where, for a comparison 
with several previously published papers on spectral ages of 
other radio sources, we have assumed the proton to elec- 
tron energy ratio to be unity. Therefore, if both values of 
_B q(class) and B oq (rev) are less than Bic/V3 (or r < r*; 
see Fig. 11), then the ages derived with the classical formula 
(open circles in Figs. 2-9) are lower than those derived with 
the revised formula (full circles in Figs. 2-9). On the con- 
trary, if both values of the magnetic field are either close to, 
or greater than Bic/V3, the 'classical' ages are greater than 
the 'revised' ages. 

4.2 Injection spectral indices 

The injection spectral indices which have been estimated 
from fits to the spectra of the lobes or the entire source 
using our measurements as well as low-frequency flux 
density values at 151 MHz from the Cambridge sur- 
veys (Green 2002 and references therein) and the VLA 
Low- frequency Sky Survey (VLSS; Cohen et al. 2007; 
http://lwa.nrl.navy.mil/VLSS) estimates for the stronger 



sources. o?i n j varies from ~0.55 to 0.88 with a median value 
of ~0.6 for our sources. For strong, non-relativistic shock 
in a Newtonian fluid cvinj = 0.5 (Bell 1978a,b; Blandford & 
Ostriker 1978). However, for relativistic shocks or shocks in 
which fields and relativistic particles alter the shock dynam- 
ics, values of ai n j are in the range of 0.35 to 0.65 (Heavens 
1989; Kirk & Schneider 1987; Drury & Volk 1981; Axford, 
Leer & McKenzie 1982). Carilli et al. (1991) estimate a inj 
to be 0.5 for Cygnus A while Meisenheimer et al. (1989) ob- 
tain similar values for a few powerful radio galaxies. For the 
sample of sources studied by Liu et al. (1992), the values 
of CYinj have been estimated from the low-frequency spec- 
tral indices using data between 38 MHz and 1 GHz. Their 
values of a ln j range from 0.65 to 1 with a median value of 
0.77. The values of the spectral indices of the hotspots es- 
timated from low-frequency measurements by Leahy et al. 
(1989), which is similar to ai n j, ranges from 0.64 to 1.16 
with a median value of ~0.8. Katz-Stone & Rudnick (1997) 
also estimate a steep injection spectral index of 0.8 for the 
compact steep spectrum source 3C190. It is important to de- 
termine the spectra of the lobes over a wide low-frequency 
range using the GMRT and upcoming instruments such as 
the Long Wavelength Array (LWA) and the Low Frequency 
Array (LOFAR) to get better estimates of the injection spec- 
tral indices, and also compare these with estimates from 
the inverse-Compton scattered x-ray spectrum (e.g. Harris 
2005). 

We have investigated any possible dependence of «i n j 
on source size using our estimates and those of Liu et al. 
(1992) and Leahy et al. (1989). We find a inj to be strongly 
correlated with luminosity (Fig. 12, upper panel), consistent 
with the well-known spectral-index luminosity correlation 
for extragalactic radio sources (e.g. Laing & Peacock 1980, 
and references therein). However, lower frequency measure- 
ments than those we have used using instruments such as 
the GMRT, LOFAR and LWA would be extremely useful to 
better determine the values of Qi n j. Omj also appears to be 
strongly correlated with redshift (Fig. 12, lower panel), but 
since most of our sources come from a reasonably narrow 
range of flux density it is difficult to distinguish whether 
the fundamental correlation is with luminosity or redshift. 
It is worth noting that while exploring correlations of ai n j 
with either luminosity or redshift, K-correction factors due 
to cosmological redshifts should not play a significant role. 
A plausible scenario for such a correlation could be multi- 
ple shocks in the heads of the jets as suggested by Gopal- 
Krishna & Wiita (1990). We have also explored possible cor- 
relation of cvi n j with source size and find that our estimates 
for the GRSs are smaller than the corresponding values for 
the smaller sources from the Liu et al. and Leahy et al. sam- 
ples, which tend to have more prominent hotspots (Fig. 13, 
upper panel). The inverse correlation of luminosity and size 
for these samples (Fig. 13, lower panel) which define the 
approximate upper envelop of the luminosity-size diagram 
(e.g. Kaiser, Dennett-Thorpe & Alexander 1997; Blundell, 
Rawlings & Willott 1999; Ishwara-Chandra & Saikia 1999) 
are consistent with the above trends. 



5 CONCLUDING REMARKS 

We summarise briefly the main results of our analysis. 
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Figure 12. The injection spectral index, Oi n j, as a function of 
radio luminosity at 1.4 GHz in units of W Hz -1 (upper panel) 
and redshift (lower panel). 
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Figure 13. The injection spectral index, «i n j, as a function of the 
largest linear size in kpc (upper panel) and the luminosity— linear 
size diagram for the sample of sources (lower panel). 



(i) The maximum spectral ages estimated for the de- 
tected radio emission in the lobes of our sources range from 
~6 to 36 Myr with a median value of ~20 Myr using the clas- 
sical equipartition fields. Using the magnetic field estimates 
from the BK formalism the spectral ages range from ~5 to 
38 Myr with a median value of ~22 Myr. These ages are 
significantly older than those of smaller sources (e.g. Leahy 
et al. 1989; Liu et al. 1992). The GRSs are older sources 
possibly evolving in underdense regions of the intergalac- 
tic medium (cf. Cotter 1998; Mack et al. 1998; Machalski, 
Chyzy & Jamrozy 2004), although the external medium is 
often asymmetric on scales of ~1 Mpc (e.g. Ishwara-Chandra 
& Saikia 1999; Konar et al. 2004, 2007). 

The spectral ages depend on the estimated magnetic 
fields. Strengths of the equipartition field calculated with 
the revised formalism proposed by BK are greater than the 
corresponding values provided by the classical formula of 
Miley (1980) by a factor of ~3 (Paper I). However, the in- 
ferred synchrotron age is a function of the ratio B jB\c- If B 
is greater than Bic/V3, the 'revised' age is lower than the 
'classical' age. 

(ii) The injection spectral indices range from 0.55 
to 0.88 with a median value of ~0.6. Our estimates for 
the GRSs are marginally smaller than those estimated for 



smaller sources by Leahy et al. (1989) and Liu et al. (1992). 
Reliable low-frequency measurements of the lobes using in- 
struments such as GMRT, LWA and LOFAR are required 
to get more reliable estimates of the injection spectral in- 
dices. We have explored possible correlations of a ln j with 
other physical parameters and find that it appears to in- 
crease with luminosity and redshift, but shows an inverse 
correlation with linear size. 
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